cholesky_rank1_downdate

public interface cholesky_rank1_downdate

Module Procedures

private subroutine cholesky_rank1_downdate_dbl(r, u, work, err)

Computes the rank 1 downdate to a Cholesky factored matrix $A = R^T R$ such that $A_1 = A - \vec{u} \vec{u}^T$ . This operation only works if the new matrix $A_1$ is positive definite.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=real64),	intent(inout),		dimension(:,:)	::	r	On input, the N-by-N upper triangular matrix $R$ . On output, the updated matrix $R_1$ .
real(kind=real64),	intent(inout),		dimension(:)	::	u	On input, the N-element vector $\vec{u}$ . On output, the rotation sines used to transform $R$ to $R_1$ .
real(kind=real64),	intent(out),	optional,	target, dimension(:)	::	work	An optional argument that if supplied prevents local memory allocation. If provided, the array must have at least N elements. Additionally, this workspace array is used to contain the rotation cosines used to transform $R$ to $R_1$ .
class(errors),	intent(inout),	optional,	target	::	err	An error object to report any errors that occur.

private subroutine cholesky_rank1_downdate_cmplx(r, u, work, err)

Computes the rank 1 downdate to a Cholesky factored matrix $A = R^H R$ such that $A_1 = A - \vec{u} \vec{u}^H$ . This operation only works if the new matrix $A_1$ is positive definite.

Arguments

Type	Intent	Optional	Attributes		Name
complex(kind=real64),	intent(inout),		dimension(:,:)	::	r	On input, the N-by-N upper triangular matrix $R$ . On output, the updated matrix $R_1$ .
complex(kind=real64),	intent(inout),		dimension(:)	::	u	On input, the N-element vector $\vec{u}$ . On output, the rotation sines used to transform $R$ to $R_1$ .
real(kind=real64),	intent(out),	optional,	target, dimension(:)	::	work	An optional argument that if supplied prevents local memory allocation. If provided, the array must have at least N elements. Additionally, this workspace array is used to contain the rotation cosines used to transform $R$ to $R_1$ .
class(errors),	intent(inout),	optional,	target	::	err	An error object to report any errors that occur.

cholesky_rank1_downdate Interface

Contents

Module Procedures

public interface cholesky_rank1_downdate

Module Procedures

private subroutine cholesky_rank1_downdate_dbl(r, u, work, err)

Arguments

private subroutine cholesky_rank1_downdate_cmplx(r, u, work, err)

Arguments